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In this article we review recent advances in the area of QCD factorization. We begin with a brief outline of 
QCD factorization, as well as a description of some recent results regarding B — > VV decays. We examine what 
is necessary at NLO, look at recent advances towards that goal and consider what remains to be done. 



1. Introduction 

Charmless hadronic two-body B decays are a rich 
source of flavour physics information and are particu- 
larly well suited to the study of CP violation. The 
theoretical description of these decays is hampered 
by the need to separate the long distance physics 
of hadronization from the short distance electroweak 
physics of interest. QCD factorization (QCDF) 0,0, 
[3] is a theoretical framework suitable for the study 
of these decays. At leading order (LO) in l/m b , it 
allows for a clean separation of scales; the long dis- 
tance physics is parametrized in terms of form factors 
and decay constants which can be related to other 
processes, while the short distance physics, including 
strong phases, is perturbatively calculable. The effec- 
tive field theory formalism necessary to go beyond LO 
in l/m& is the Soft Collinear Effective Theory (SCET) 
developed by Bauer, Fleming, Pirjol and Stewart. [4[ 

QCDF was supplemented by an effective field the- 
ory description, SCET, giving further credibility to 
the factorization result first put forward in the con- 
text of QCDF. A modern view of QCDF, that we 
adopt in this short review on the status of QCDF, 
would be that QCDF is the application of SCET to 
two-body B decays. We therefore use the terms QCDF 
and SCET interchangeably throughout the review. In 
the application of the theory, two approaches exist: 
the original approach of Bcnckc, Buchalla, Neubert 
and Sachrajda (BBNS), often called QCD factoriza- 
tion (QCDF) and the approach advocated by Bauer, 
Pirjol, Rothstein and Stewart (BPRS) often referred 
to as SCET. Although this article will make reference 
to both approaches, emphasis will be placed on that 
of BBNS. The discussion reflecting recent results from 
the BPRS perspective has been covered at this confer- 
ence by Christian Bauer in his presentation on SCET. 

In the next section we give a basic description of 
charmless two-body hadronic B decays in the QCD 
factorization formalism and describe the differences 
which exist between the BBNS and BPRS approaches. 
In Section [3 we discuss a recent application of QCD 
factorization to B — > VV decays. In Section^] we dis- 
cuss the theoretical situation at NLO, including both 
l/mb and Oi s corrections, some recent work towards 
the NLO perturbative corrections and what work re- 
mains to be done. 



Figure 1: A typical charming penguin contribution 
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2. The Basics of Factorization 

At LO in l/m b , the amplitude for B — > MiM 2 de- 
cay factorizes, 

(M 1 M 2 \0\B) = F B ^ Kh T l * /m 2 $m 2 + Mi «-» M 2 

+T U * f B §B * fM^Ah * fM 2 ^M 2 - (1) 

T l and T n are hard scattering kernels that depend 
on the specific process being considered, but are 
perturbatively calculable in an expansion in ot s (mb) 
and a s (fii) , where [ii is the intermediate scale fii = 
i/AqcdW6- Jai and fg are light and heavy decay 
constants, while <I>m and $b are light cone distribu- 
tion amplitudes (LCD A). They are universal and do 
not depend on the decay process. They are deter- 
mined from QCD sum rules, lattice QCD, or exper- 
iment. Although this formula is valid only at LO in 
1/mi,, it is valid to all orders in a s . A systematic 
treatment of the higher order corrections can be given 
using SCET. 

2.1. The BBNS and BPRS Approaches 

As discussed in the introduction, there are two main 
approaches which we refer to as BBNS P, 0, El and 
BPRS Below, we outline several of the important 
differences, which have been discussed at great length 
in HQ. 

Perhaps the most important difference between 
the two approaches is the treatment of the so-called 
charming penguin contributions (see Fig. fJJ. BPRS 
argue that since 2m c ~ mb, there are configurations 
with almost on-shell charm quarks leading to long dis- 
tance effects. As a result, they parametrize the pen- 
guins and fit them from data. BBNS maintains that 
charming penguin diagrams factorize in the formal 
l/mb limit and should not be given special treatment. 
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In the BBNS approach a subset of higher or- 
der 1/m.b corrections are typically included that are 
deemed to be numerically large. This is in contrast to 
BPRS, where calculations are consistently performed 
to a fixed order in the 1/mb expansion. 

It has been the philosophy of BBNS that the inputs 
used to make predictions for B — > M1M2 decays do 
not come from these decays themselves; rather, the 
form factors and decay constants enter from other ob- 
servables, lattice calculations, or some other external 
source. BPRS on the other hand impose no such re- 
striction and use not only external sources but also the 
B — > M1M2 decays themselves as inputs. In particu- 
lar, BPRS has up to this point chosen to fit the SCET 
form factor, S (see below), rather than perturbatively 
expanding it and relating it to other observables. 

BBNS and BPRS also differ in the counting of the 
soft overlap and hard scattering terms in (JTJ. In 
both cases these terms are counted as leading order 
in l/m h , but BBNS assumes that soft overlap is dom- 
inated by soft contributions without a s (/Zi) suppres- 
sion, and is thus assumed to be larger. 

The differences given above leave the BBNS ap- 
proach more predictive than BPRS, but at the cost of 
being less conservative. We will see that when NLO 
corrections are taken into account, some of these dif- 
ferences will disappear. 



3. Recent Results for B — > VV 

QCDF had been used to predict the branching ra- 
tios and CP asymmetries of a large number of two- 
body charmless B — > PP and B — > PV decay modes 
In a recent work, Beneke, Rohrer and Yang 
have examined B — * VV decays in the context of 
QCD factorization. 

According to naive factorization, the ratio of longi- 
tudinal, negative and positive helicity B — > VV decay 
amplitudes scale as [l(| 



A 



A :A_:A + = 1 : — :[ — ) 



A V 



(2) 



Although this scaling is well obeyed by B — > pp de- 
cays, it is badly broken by penguin dominated decays 
such as B — + K*p and B — ► K*<fi, which have anoma- 
lously large transverse polarizations [Til ] 

Beneke, Rohrer and Yang have identified a novel 
contribution to B — ► VV decays that had not been 
previously taken into account. This contribution, a 
contribution to the transverse electromagnetic pen- 
guin amplitudes, is given by the diagrams in Fig. [21 

Although this contribution is in fact 0(m b ) en- 
hanced over the leading terms, it is also a cm sup- 
pressed. As a result, it is the dominant contribution to 
the transverse helicity modes, but is not large enough 



Figure 2: New contributions to B — » VV decay 
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to explain the puzzle of the large transverse polariza- 
tions of K*p and K*<f>. 

Although insufficient to resolve the discrepancy be- 
tween theory and experiment, the new contribution 
has pronounced effects, as demonstrated in the exam- 
ples below. These examples make use of the current 
experimental data to estimate the QCD penguin con- 
tributions to the decay amplitudes. QCD factoriza- 
tion is used to estimate all other contributions, includ- 
ing the electroweak penguin effects. (Several other ap- 
proximations are also necessary as explained in 0.) 

As a first example, consider the ratio of CP aver- 
aged helicity decay rates 



js*o\ 



Sh 



2T h (p a K 
f h (p-K*o) 



1 



pEW 



Ph 



A h , (3) 



where h denotes the polarization (h = 0, — , +), P^ w 
is the electroweak penguin contribution calculated in 
QCDF and Ph is the QCD penguin contribution whose 
magnitude is determined from the Br and /l of B — ► 
p~K*° (note that A h (p~ K*°) = P h ). A h is a small 
contribution from the color-suppressed tree amplitude 
(with some dependence on P h EW ). Without the effects 
of Fig. |3 SL takes on the value 0.7 ± 0.1, but is 
increased to 5_ = 1.5 ± 0.2 when the new effects are 
introduced. 

Also consider the ratio 



_ 2T h (p°K*-) 
h f h {p-K**) 



pEW 



Ph 



(4) 



As before, \Ph\ is determined from the Br and fi, for 
B -> p-R*°. This analysis leads to 5"- = 0.5 ± 0.1 
with the new contribution and S- = 1.2±0.1 without. 

Here it is possible to compare this result to ex- 
periment. The above ratio of branching ratios is 
related to the ratio of polarization fractions f' h = 
fh(p°K*-)/fh(p-K*°) so long as the CP asymme- 
try is assumed to be zero. Using theory to estimate 
this asymmetry leads to the results 



Jo = 1.3 ±0.1 
f_ = 0.4 ±0.1 



[1.1 ±0.1 (without)] 
]0.8±0.1(without)]. 



(5) 



The experimental results for these ratios are JT2I IH 



ULUdLdlia fo = 1.45±8;£ and f_ = 0.12+"^. As 
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can be seen, the effects of the new operator can be 
large (as is the case for f'_) and would seem to bring 
theory further in line with experimental expectations. 



4. NLO Corrections in a s and l/m b 

In this section, we examine current status and out- 
look for the NLO corrections to the factorization for- 
malism (be it SCET or QCD factorization). The NLO 
corrections include both perturbative corrections to 
the calculable kernels (T l and T 11 ), as well as power 
corrections {l/mb) to the leading order result Eq. (fj^). 

The primary motivation to include the NLO correc- 
tions is simple: to reduce overall uncertainties. Both 
perturbative and power corrections are expected to 
be of order ~ 20%. Thus, to reduce theory uncertain- 
ties below this level and remain competitive with ex- 
periment, it becomes necessary to include both these 
corrections. 

There are also observables for which the NLO cor- 
rections are indispensable. CP asymmetries will re- 
ceive large corrections at NLO, if at LO the "tree" am- 
plitudes (for definition see e.g. 0) do not originate 
from the Q\ t i operators of the effective electroweak 



Hamiltonian. For AS = 1 modes such as B® 



<>„(') 



tt r] 



or B — > K°ir we expect large 0(1) corrections to 
their CP asymmetries. This situation is even more 



extreme for AS = modes such as B 



K V K- or 



B° — > K°K°, where no CKM hierarchy exists between 
the "tree" and "penguin" amplitudes. 

4.1. Status of NLO corrections 

Typical power corrections are of the order ~ 20%, 
and will likely introduce many new parameters. SCET 
is the theoretical framework necessary to describe the 
l/mb corrections, but almost no l/mb corrections 
have been considered. Exceptions include the anni- 
hilation contributions and chirally enhanced contri- 
butions that have been taken into account by BBNS. 
Much work remains to be done for these corrections. 

Unlike the power corrections above, almost all of the 
order a s corrections have been calculated. To exam- 
ine these corrections, it is useful to rewrite the basic 
QCDF formula (TTJ, by writing T H = H H * J, 

(MiM 2 \0\B) = ( F B^M 1T I + E B^Ah ^jjllj 

*/a/ 3 *m 2 +Afi <-> M 2 (6) 



where 



—B^Mi 



J * fs^B * 



(7) 



T l and T u are hard coefficients, calculable in 
a s (mi) ~ 0.2 where perturbation theory is expected 
to be valid. The LO contribution to the hard coeffi- 
cients is in general 0(1). 



The jet function, J, is perturbatively calculable in 
ot s (jii) and differs from the hard coefficients in two im- 
portant ways: the LO contribution to J is 0(a s (//j)) 
and J depends on only one of the two final state 
mesons, and so can be related to other processes. 
In particular, the same jet function occurs in semi- 
leptonic decay. As of the last FPCP conference, FPCP 
2004, the hard kernel T 1 was known to 0(a;?(m&)/3o) 
(NLO+) [H |H |2j , H u was known to (a° (m 6 ) ) 
(LO), and J was known to O(o%(m)) (NLO). 

4.2. Corrections to H 11 

Recently a subset of the 0(a s (mbY) (NLO) cor- 
rections to H u have been calculated [23. This cal- 
culation involved matching the effective electroweak 
Hamiltonian onto leading order SCET operators at 
one loop. In this calculation, only the operators Qi,2 
of the electroweak Hamiltonian were matched, leav- 
ing the remaining operators for a future work. Taken 
with the earlier results for the NLO calculation of the 
jet function J, « s ^ M i became known to NLO for a 
subset of operators. 

As a phenomenological application of these new cor- 
rections, the authors re-calculated the predicted tttt 
branching ratios. The amplitudes for these decays are 
given by (C = i^m%f n f?*(0)) 



y/2A B —*-ifl = c[v ub v: d 



Q'l + Q'2 



a 1 + Q'4 



Q'2 - Q4 



+ v cb v; d &i}, 

V cb V: d a%), 



where the ai are the amplitude coefficients defined in 
Q. Only the tree amplitudes 01,2 are affected by the 
new contributions; they become, 



a% {tttt) 



1.015 + [0.025 + 0M2i] v - [0.009] tw3 
-0.014lo- [0.024 + 0.020i] N LO 



n qqo+0 029 



-O.OOTjJSS)*, 



(8) 



0.184 - [0.152 + 0.077i] v + [0.056] tw3 
+ 0.088 LO + [0.029 + 0.034i] NLO 



o.205ig:iS + (-o.o43i8:gg)i. 



(9) 



Here the leading order (tree) terms are given no sub- 
script, the NLO (0(a s )) vertex corrections are de- 
noted with a subscript V, the included subset of power 
corrections are denoted with the subscript tw3, while 
LO and NLO denote the LO and NLO (new) spec- 
tator scattering corrections. As can be seen, the new 
corrections contribute important new strong phases to 
the amplitudes. 

The above numbers offer an important insight into 
the validity of perturbation theory at the intermediate 
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Table I CP-averaged branching fractions (xlO -6 ) Val. 
denotes the central theory prediction, CKM, hadr. and 
pow. denote the theory uncertainty associated with 
CKM matrix elements, hadronic parameters and power 
corrections respectively, and expt. denotes the present 
experimental value. 



Mode 


Val. CKM hadr. pow. 


expt. 


B~ -> 7T~7r° 
B° -> 7V + 7T- 

B° -» 7r°7r° 


r r +0.3 +0.5 +0.9 
°-° -0.3 -0.4 -0.8 
c n +0-8 +0.3 +1.0 
°- u -0.9 -0.5 -0.5 
n 7 +0.3 +0.5 +0.4 
u -' -0.2 -0.2 -0.3 


5.5 ±0.6 
5.0 ±0.4 
1.5 ±0.3 



scale i/AQCDWf,. In ev2 (""""), the NLO spectator scat- 
tering terms are smaller than the LO terms in agree- 
ment with naive expectations. In ai(7T7r), the NLO 
corrections appear to be enhanced over the LO term, 
and may lead one to believe that perturbation theory 
is breaking down. This is not a valid conclusion how- 
ever as the LO term is suppressed as it is proportional 
to a small Wilson coefficient. These corrections there- 
fore suggest that perturbation theory is valid at the 
intermediate scale. 

In trying to reproduce the 7T7t branching ratios (see 
Table U), the data favours inputs where 



Btv 



\v ub \n 



(0) 



8.10 x 1CT 4 - 0.775 [\V ub \f^( 

/tt/b 



m b f^(0)X B 



= 1.96 



m b /f 7r (0)A i 



I def ' 



(10) 



def 



with the subscript def. refering to the default values 
used in the determination of ai^mr) above. These 
inputs are consistent with values of the form factor 
and Xb smaller than their default values 23], 



/£*(()) = 0.22 [0.28±0.05] def , 
X B = 0.23 [0.35±0.15] dc f, 



(11) 



and brings the form factor closer to the BPRS value: 
0.176 ± 0.036 0, [H|. As well, these values remain 
consistent with what is found using lattice and semi- 
leptonic data p3 |. 



4.3. At NLO, BBNS - BPRS = ? 

It is instructive to ask what differences remain be- 
tween the BBNS and BPRS approachs once NLO cor- 
rections are taken into account? At least one of the 
differences is no longer present at NLO. In the BPRS 



was 



approach at LO, the SCET form factor E B ^ Ml 
fit from B — ► M\M^ data. At NLO one can no longer 
fit for j there is a non-trivial convolution and 

not simple moments as at LO. One is forced to follow 
the path of BBNS and expand £ B ^ M i in terms of the 
jet function J as well as other long distance factors. 

The question as to whether the soft overlap func- 
tion should be fit from B -> MM data (BPRS) or 
be determined by use QCD sum rules or lattice QCD 



(BBNS) remains unresolved. BBNS may be moving 
in the direction of BPRS in that the value of A seems 
to be determined from fit to B — > 7777 modes. 

One issue that remains unchanged is what to do 
with the charming penguins. For the near future, the 
disagreement as to whether the charming penguins 
should be considered to have long distance effects and 
fit (BPRS) or whether they are assumed perturbative 
and calculated (BBNS) will remain an open question. 
Fitting the charming penguins is always a valid pro- 
cedure, but lacks the predictive power of perturbative 
calculation. A comparison of theoretical predictions 
between the two approaches should one-day resolve 
the issue. 



5. Conclusions 

QCD factorization is an important tool in the study 
of flavour physics. Reduction of the theoretical un- 
certainties associated with QCD factorization predic- 
tions requires NLO corrections to be computed. Al- 
though the complete set of NLO in a s corrections is 
almost complete, only a handful of l/m?, corrections 
have been considered. 

In this work we presented a brief summary of the 
QCD factorization formalism as well as some recent 
developments and discussed its similarity to SCET. 
We examined several of the differences between the 
approach put forth by Beneke, Buchalla, Neubert and 
Sachrajda (BBNS) and the approach advocated by 
Bauer, Pirjol, Rothstein and Stewart (BPRS). We 
concluded that many of the characteristics that differ- 
entiate between the two approaches are not applicable 
when NLO corrections are taken into account. 
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